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Abstract. Motivated by the surface of topological insulators, the Dirac anomaly's 
discontinuous dependence on sign of the mass, ra/|ra|, is investigated on closed 
topologies when mass terms are weak or only partially cover the surface. It is 
found that, unlike the massive Dirac theory on an infinite plane, there is a smoothly 
decreasing current when the mass region is not infinite; also, a massive finite region 
fails to exhibit a Hall current edgeexerting an extended proximity effect, which can, 
however, be uniformly smalland oppositely orientated Hall phases are fully quantized 
while accompanied by diffuse chiral modes. Examples are computed using Dirac energy 
eigenstates on a flat torus (genus one topology) and closed cap cylinder (genus zero 
topology) for various mass-term geometries. Finally, from the resulting the properties 
of the surface spectra, a potential application for a flux-charge qubit is presented. 
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1. Introduction 

The 2+1 dimensional massive Dirac theory on an infinite plane exhibits an anomalous 
contribution to the Hall conductance (h, c, = 1, e = |e|) cr xy = jj^f^ [UEIE] representing 
a /ia// charge pumped contribution per unit of threaded flux in a Laughlin type 
setup [U |5]. Unfortunately, this half quantization is not directly observable, since a 
at two-dimensional lattice theory such as that for graphene, must have an even number 
of Dirac modes [6J, [TJ, [8] . Having effectively two or more pairs of fermion flavors, single 
fermion results are effectively doubled. 

In contrast, the surface of 3-dimensional topological insulators contain an odd 
number of massless Dirac fermion modes on closed topologies [9j dUJ [EU [[21 CESJ [HI 
[T5l [T6l [T71 [EE]. Effective Dirac mass terms could be given to the surface theory at 
selective regions, by coupling to a ferromagnetic time-reversal breaking material [19J. An 
application of the naive infinite plane results to this system suggests several questions. 
If the massive region only covers a small area, will charge accumulate on its edge? If not, 
can a sign change of a weak local mass induce a sign flipping transition of the quantized 
current globally? Finally, how do magnetic fluxes induce fractional charge pumping in 
specific cases? In the present work these issues are investigated and clarified. 

Rosenberg et al [20] gave a solution to the last question by concluding that fractional 
charge accumulation did not occur when a flux is inserted through a topological 
insulator. The surface theory was said to breakdown at n flux, as gapless bands induced 
in the bulk were responsible for discharging the ends of the flux tube. This was called 
the wormhole effect. Yet, it is not possible to explain a half-charge quantized transport 
on the basis of bands. In fact, the two dimensional plane result implies the existence 
of the current in a fully gapped theory. The results to be presented will clarify these 
points and it is conjectured that fractional charge could indeed be observed. 

In the present work it is emphasized that the surface theory is sufficient to resolve all 
previous questions. Three results will be shown through solved examples and motivated 
from the effective action of topological insulators. These are: 1) that mass terms 
do not induce perfect quantization in weak limits for closed finite surfaces (partial 
quantization), 2) that mass edges do not result in fractional Hall conductance edges 
("extended proximity" effect) and 3) when the transition region separates oppositely 
orientated masses those phases are fully quantized but overlapped by a sufficiently diffuse 
chiral band. The existence of chiral bands is not a new result [TSl ED E2], but rather 
the full quantization of the separated phases. It is also suggested that the wormhole 
effect [20] only occurs in the localized flux limit and in that case can be thought of as 
an example of the extended proximity effect on a genus one surface. 

Section [2] will first offer motivation for these results from the effective surface action. 
Sees. [3] and [4] will then present specific examples and constructions showing how the 
effects are manifested by the surface Dirac theory. The solutions to the Dirac equation 
has been investigated in other contexts and geometries involving topological insulators 
and graphene [18j [23J EH E5]. The anomalous Hall current is understood from the 
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wavefunction spectral asymmetry (see Appendix A). Finally, gaining intuition from the 



different mass geometries computed, a potential configuration of masses with inserted 
flux is proposed as an architecture for a qubit in section [5j 

2. Effective Action 

The conclusions reached in this work can be motivated from the effective field theory 
for topological insulators. Starting from the topological BF theory it was shown in 
reference [26] that the entire electromagnetic response comes from the surface. After 
integrating out the effective fields, the Hall action is obtained: 

Ssurf = ~ J e^eA,d v A p . (1) 

Here, and throughout, Greek indices take values over the 2+1 space-time describing the 
surface. 9 is a parameter or background field giving a hall conductance a xy = #e 2 /47r 2 . 
9 is required to be ±7r to respect time-reversal invariance, and comparing with the Dirac 
theory on a plane allows one to identify 9 =Sign(m)7r [TJ [2j [3]. 

Partial Quantization. First it is noted that on the surface of the insulator there 
is no reason to expect 9 to be quantized if time-reversal invariance is being broken by 
mass-terms. Importantly, for the Abelian theory on a closed surface (as in equation [I]) , 
gauge invariance is consistent with 9 unquantized. A fixed value of 9 other than ±7r 
would lead to a an unquantized Hall current. 

Extended Proximity Effect. For a more complicated arrangement of mass terms it 
might be allowed 9 — » 9{x) ) as an effective parameter not necessarily equivalent to the 
local sign of the mass. An abrupt change in the value of 9 from non-zero to zero can be 
thought of as a Hall edge. However, gauge invariance in general requires that 9 remain 
constant: 

S surf - J e^dd^kdvAp = j €*">Adp6d v A p ^ 0, (2) 

(0 mod 2tt implied for the last term). Therefore 9 ^ 9{x) which is presumably 
responsible for the extended proximity effect and the absence of a Hall edge. 

Quantization of opposite phases. The exception to equation [2] is if a region exists 
of zero magnetic field (which is the case for closed surfaces that do not surround 
monopoles). In these cases 9 can transition (over a region of zero field) from one fully 
quantized (half-charge) Hall phase to the oppositely quantized phase accompanied by a 
chiral band. 

To see this, working in the Coulomb gauge (A = 0), zero field strength requires 
Ai = diT (i G 1,2). Inserting this into the action where 9 varies, naively gives: 

o2 



Se&o = —, I e^d^doTdjT). (3) 



e 

47T 2 

Normally the electromagnetic field is physically determined by a single choice of T. Then 
there are no unique modes and no way to form propagating packets. However on the 
compact surface there exist different disjoint T distinguished by winding number. These 
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different winding numbers can act as different modes forming a compact (discrete) chiral 
band. 

If 9 varies in the j direction, the i-direction is topologically equivalent to a circle 
separating the two phases, and with periodic boundary conditions in time, the action 
can be integrated to give (setting e = 1): 

S 9 # = A0kl (4) 

where fc, /, are integers. This remains gauge invariant for A9 = 2nn (n G Z). Therefore 
9 must separate two fully quantized oppositely orientated phases. The complete picture 
is that the chiral band can absorb or release a whole charge, half of which comes from 
each separate phase. The modes may be diffuse over the entire region where the field 
strength vanishes, however. 

In the rest of this work, these suggested results will be shown in various solvable 
examples. 



3. Flat Torus 



The first easily solvable case is a flat torus, which can be considered as a simplification 
for a crystal wafer with a hole drilled through it. Labeling coordinates z and (j) for the 
two orthogonal directions on the torus, a flat metric can be used. It is important to note 
that four inequivalent forms exist for the Dirac equation on the torus, corresponding to 
the inequivalent spin structures. For example writing the torus as a square with edges 
identified naturally leads one to a Cartesian-like form (with constant Pauli matrices). 
However, this will not agree from an embedding in 3-dimensional Euclidian space as used 
in references [27l l28l [20] . While both forms are mathematically consistent, they lead to 
different physical results. The spectrum, for example, has a finite-size gap present in 
one case but not in the other. However, none of the conclusions made here will depend 
on this choice. Therefore I will pick the spin structure that reproduces the half-integer 
azimuthal dependence found for the unique choice in the case of a closed (with caps) 



cylinder (see section [4] and Appendix B). It is also the choice which leads to a derivation 
of the Dirac equation by embedding. 



3.1. Weak Hall current 

To ascertain an anomalous Hall conductance as a function of mass-strength I consider 
a fully massive Dirac fermion with a step potential in the (j) direction V((j>) = 
v(Q((j> + 1/2) — 0(0 — 1/2)) — v/2, corresponding to localized electric fields at ±//2, 
E = [—S((j) + 1/2) + 5{(f) — l/2)]v(j>. This construction is shown in figure [lj The units 
are set with L\ = 2n and other quantities are measured in units of Li, L 2 —> 27rL 2 /Li, 
I -> 2nl/L\. Then the Dirac equation is (7 = a 3 ,Y = aV 1 ^ = a 3 a 2 ): 

e- % ^{-io Y d z - ia 2 d^ + a 3 m - eV^e^^ = E^j. (5) 
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Figure 1. (Color online) Coordinates and dimensions (left) for the smooth torus 
(right). Note sides are identified in the left figure. The surface is massive (constant) 
everywhere, but a non-zero scalar step-like potential is introduced with values as shown 
by V (right), z is in the vertical direction at right. Note the torus is separated into 
two regions, or three with edges identified, which will be labeled as I, II, III. 




Figure 2. (Color online) Spectrum for the configuration of figure [T] with m = 10, 
v = 4, / = 2, ^2=5. Units are in terms of 2ir / L\ as defined in the text. 



with < z < 2tt and — L 2 /2 < cj) < L 2 /2 while mass and energy are normalized to units 
of 2n/Li. The solutions in each region are matched and the quantization condition 
derived in |Appendix B[ The corresponding spectrum is unremarkable and shown in 
figure [2] as a function of the azimuthal angular momentum. It remains fully gapped for 
any threaded flux. 

Nevertheless the anomalous Hall current, j z {(j)) ) is non-zero and can be computed 
numerically using equation ( |A.5[ ). According to the infinite plane result, one expects 
3 —^- = [5((f) + l ( f ) /2) — 8(<j) — /^/2)]|Sign(m). However, the computed result turns out to 
have a smoother mass dependence as shown in figure [3| illustrating partial quantization. 
The current smoothly approaches zero as the mass approaches zero with the onset 
of full quantization only with a sufficiently large mass-distance m ~ 20 (in units of 
2n/Li). If the current is computed along any window other than the delta function 
peaks, none is found, as expected. There is no further dependence on the potential 
strength beyond the expected (linear) scaling. Finally, note that the system exhibits a 
current despite remaining fully gapped; consistent with the understanding that bands 
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Figure 3. (Color online) Hall current, / (in the ^-direction) in units of — e/2v, for the 
configuration of figure [l] along a strip of width = L2/IOO centered around (j) = —1/2. 
From the infinite plane result, 1 = 1 expected. The x-axis labels the cutoff energy 
employed (that is, all states below said energy are summed to approximate the delta 
function peak). Each line represents a different mass: from top to bottom m = 40, 20, 
10, 5, 2.5, .5 (in units of 2i\ /L\). All cases take v = 4, I = 2, and L2 = 5. The same 
values with opposite sign are obtained along the strip <\> = 1/2. 



are not responsible for fractional charge transport. 
3.2. Extended proximity effect: first suggestions 

A second case that can be easily computed suggests that while the current depends on 
the mass strength, it does not respond locally to a massive region. To see this, a torus 
with a strip of mass m is considered as shown in figure [4j It is checked that no charge 
accumulates after flux threading, in other words, a strip of mass does not behave like an 
edge and a Hall phase is globally induced. While only suggestive, much more compelling 
examples will be considered in the next section. 

The Dirac equation and solutions in this case are given in |Appendix B and 



eigenvalues are shown in figure [5j The large (blue) points are at the allowed azimuthal 

quantum numbers k G Z + 1/2 of the wavefunction ip defined in the appendix as 

~ 3 

ip a = e l ^^ilj a . Because of the half-integer quantization there is a finite-size energy 
gap, of order 1/L 2 . When a flux <£> is inserted its effect can be undone through a 
transformation which amounts to shifting k k + $. The smaller (red) points show the 
evolution of the eigenvalues under this threaded flux. 

For the massless case (figure [5](a)), the bands traverse the gap. These bands are 
not responsible for the Hall current and half-integer charge flow around the torus. The 
wavefunction for these states satisfy ipE<o = & 3 ^e>o =constant which does not support 
a current in the z direction or charge accumulation. There is a constant current in the 
(j) direction = 27re 2 , which is the usual longitudinal current expected of free fermions 



Anomalous Currents on Closed Surfaces 



7 




Figure 4. (Color online) Coordinates and dimensions (left) for the smooth torus 
(right). Note sides are identified in the left figure and (j) and z have been exchanged 
relative to figure [I] However, at right, z is still in the vertical direction. The torus is 
separated into two regions, or three with edges identified, as shown in at left labeled 
I, II, III. A mass term, m, is present only in the shaded area, region II. 



(a) (b) 

E E 




Figure 5. (Color online) Spectrum for the configuration of figure [4j (a) massless case, 
m = with I — .5 and L 2 = 2. (b) with mass m = 1.5, and same aspect ratios. Units 
are in terms of 2ir / L\ as described in the text. 



accelerated under potential 2ne. 

For any small strip of mass, a gap forms for all fluxes as shown in figure [5](b) . In 
this case the flow returns the vacuum to itself, and again no charge accumulation is seen 
suggestive of a extended proximity effect and lack of Hall edge. This spectral flow is to 
be contrasted with the results of the next section. 

4. Surface of a Closed Cylinder 

The surface of a closed cylinder offers analytic solutions for the more interesting case 
showing fractional charge accumulation and direct detection of weak quantization and 
the absence of a Hall edge. 

4-1- Closed Cylinder Cases: weak currents and extended proximity effect 

The relevant dimensions and coordinates for the cylinder surface are described in 
figure |HJ The Dirac equation in the three different sections and matching conditions 
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Figure 6. (Color online) Closed cylinder of radius R, and side (cap-to-cap) length d. 
Three separate sections are labeled I, II, and III each with local coordinates as shown. 
The caps, I and III, have local polar coordinates with r x <j> oriented upward. <j> globally 
has the same orientation. The side has z oriented in the upward direction, with the 
origin set in the center of the side. Different cases are considered where each section 
may separately have a constant oriented mass m a . See text. 



are derived in Appendix C, they are on the top and bottom caps (with and without 
masses) : 

r 



a 1 sin < 



[(a 1 cos (j) + a 2 sin (j))d r + (a 2 cos < 

+ rai,(iii)0" 3 ]^,(iii) = #0i,(m) (r, 0) 

and on the side: 

e- i ^\-icj 1 d z - ia 2 d^ + a 3 m n )e^ fyn = Eipu(z, 
or defining in each region i/j = exp(z^-0) , 0, 

, odd. 



(6) 
(7) 



-ia l d r 



1(7 



10 
~2r 



^i,(iii)0" ^1,(111) = Eip: 



X(iii), 



[—i(j l d z — io 2 d$ + m n cr 3 ) ^ n = £^11 . 



The matching conditions are (see Appendix C): 

a 2 ^ u 



0: 



l\r=R 



\z=d/2, lfoll\r=R — Wll\z=-d/2- 

I consider six cases (m>0^ m(x))\ (a) massless case, mi = ran = ram = (b) positive 
mass on top cap alone rri\ ra, ran rani 0, (c) positive mass on the side ran — ra, 
rai ram (d) positive mass on top and bottom m\ —ram ra, (e) positive mass 
everywhere rri\ —ran —ram ra, and finally oppositely oriented masses on the two 
caps (f) rai = +ram = ra, ran = 0. Note the relative minus sign between the top and 
the rest of the cylinder in the cases (b)-(e) to describe a same sign mass, comes from the 
a 2 transformation in the top-to-side matching; stemming from the opposite orientation 
of the r and z directions at the top-to-side boundary. 



^11 1 



(8) 

(9) 
(10) 

(11) 



Anomalous Currents on Closed Surfaces 



9 



Unlike the torus, the effect of the flux is not simply to shift the azimuthal quantum 
number. In particular some boundary condition must be implemented at the origin 
which physically involves the flux details. Interestingly a simple shift alone of k 
(azimuthal quantum number of iji) would allow for extra solutions, unconstrained 
by normalizability at the origin (see reference [29]). The simplest possibility is to 
take a profile for the flux as a localized delta-function-ring with small radius, e: 
B(r) = <M(r — e)/(27re) [29J, which is expected to be qualitatively similar to a diffuse 
flux of width e. 

Figure [7] shows the spectrum for the different cases as a function of azimuthal 
quantum number k. As before, the large (blue) points represent the spectra with no 
flux, and the smaller (red) points are the evolution of the spectra as one unit of flux 
is inserted. The effect of the flux is manifested by the appearance of extra states at 
the end of the cycle inserting one unit of flux (marked A and B). The creation of 
new states relative to the original vacuum will have a half charge associated with them 
(see Appendix A). For all cases a simple pattern emerges: anomalous bands create extra 



states near —mi and mm (A and B points). One can show that the states are never 
exactly at those values, including E = for the massless case. All other bands flow 
back to, or very near, the original spectrum. The difference between the cases occurs in 
the relative distribution of the wavefunctions after one unit of flux is inserted. This is 
now described. The Fermi energy is assumed to be at zero or mid-gap. 

Partial quantization. First, in case (a) (massless case) no charge accumulation is 
seen. Indeed, the two states that appear near E — each are equally split between the 
top and bottom flux-piercing, so that no net charge density appears accumulated. In 
cases (b) through (e) , the new states are now unevenly localized resulting in a net charge 
pumped. The negative energy states are localized near the bottom cap flux while the 
positive near the top. However, the total charge pumped to(from) the cap depends on 
the mass-region strength, going to zero as this value goes to zero. Quantitative results 
for the pumped charge are shown in figure [8j showing incomplete quantization if the 
mass-region is too small. 

Extended Proximity Effect. While the positive and negative bands both approach 
the value of the local mass, the net charge pumped responds to the configuration 
globally. A striking example is case (c), where charge is pumped to the caps, despite 
the mass being only present on the side and the similarity of the spectrum with case 
(a). Also in case (b), even though the mass in present on one cap, the amount of charge 
pumped is equal and symmetric on the top and bottom caps. In case (b) the equality 
of \Qtop\ = \Qbot\ is verified to high accuracy for all the masses shown in figure [8] up to 
long cylinders of ratio d/R = 8, thereby indicating that the whole cylinder behaves in a 
similar fashion as though a mass were present everywhere. 

Opposite phases and chiral mode. An example of the expected chiral band 
separating the two Hall phases is shown in this case (f). As in the previous cases 
two new states now appear at E — — \m\. Both of these states represent a half-charge 
at each cap relative to the flux-less vacuum. After unit flux these states give a deficit 
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Figure 7. (Color online) The spectra corresponding to the cases (a) through (f) 
defined in the text. The insets reproduce cylinders for each case with color scheme: 
yellow for the massless regions, blue for positive mass, and in the case of (f), red for 
a negative orientated mass. Blue (larger) points at half-integer azimuthal number k 
are the spectra with no flux inserted. The red (smaller) points represent the spectral 
flow as a localized flux is inserted as described in the text and shown in figure [6j After 
one unit of flux, new states appear at marked points A and B. For all cases \m\R = 5, 
d/R = 1/2, and the energy, E : is in units of 1/R. 



from each cap while a chiral mode becomes occupied. 

As a final note, if the solenoid could pierce only one cap (a net outward flux), then 
a single E = ±m (with m the mass of that cap with flux) state would be found for 
$ = 1 in agreement with Index Theorems (see references within [USD]). In the present 
case the pair of new states ensure the overall system remains neutral but the relative 
distribution of charge does not remain uniform. 
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Figure 8. (Color online) The total charge x2 (Q in units of — e/2) in the top cap 
after one unit of flux is inserted; cases (b) through (e) clockwise, with (b) starting 
at top- left. The horizontal axis plots a log scale of m which is physically mR since 
units are normalized by R. The various lines are for ratios d/R — 1/32,1/8,2 for 
the circle(blue), square (purple), and diamond (yellow) respectively. The side-mass case 
also has d/R = 1 shown in triangle (green). 



4-2. The wormhole effect and surface theory fidelity 

Before concluding with the examples, I wish to comment on whether the closed cylinder 
results of the previous section are applicable for a real system. It was argued in 
reference [20] that surface electrons will tunnel through the bulk of a topological insulator 
for a very localized flux or if a hole is bored through the material. This was called the 
wormhole effect. Nevertheless, I will consider this question open for a general flux. 
Noting the results of section [3] one notices similarities with the work of reference [20j . 
While I considered a smooth case, with a finite-size gap 1/L 2 (blue points of figure [5]) 
it is clear such a gap will be of order 1/R where R is the radius of the interior of a 
bored cylinder. In these cases, the effect of inserting a flux tube through a bored hole 
or a single plaquette of the microscopic lattice, can be removed by simply shifting the 
azimuthal number and the spectral flow will be similar to figure [5](a). 

In fact, a similar picture will result whenever a simple global transformation can 
remove the effect of a threaded flux even if it goes through a hole or single plaquette. 
Thus the following ansatz is proposed: if an extremely localized flux string manages to 
pierce a single or a few plaquettes throughout the bulk, then currents can be described 
by a "surface" of genus one as in figure [91(a), since either for a lattice or continuum theory 
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(a) Microscopically localized flux 




(b) Microscopically diffuse flux 



Figure 9. (Color online) Illustration contrasting a microscopically localized flux (a), 
which is qualitatively equivalent to the genus one surface discussed in section |3j and 
a diffuse flux (b) which remains faithfully described by a surface theory of genus zero 
and penetrating flux. 



the flux simply shifts the azimuthal number. The spectrum of the surface theory of the 
torus contains the gap-closing seen in reference [20] and corresponds to the ability for 
electrons to propagate through the interior surface of the torus. As noted in section [3| 
this band is actually not responsible for the anomalous current. In any case, the fact 
that charge will not accumulate with a mass partially covering the surface, is just a case 
of the extended proximity effect that forces the Hall current not to exhibit an edge. 

If instead the flux is smoothly varying and extended over many plaquettes, then it 
is still reasonable to expect the low energy theory on the surface be described faithfully 
by the Dirac theory but on a surface of genus zero (this is illustrated in figure [9](b)). 
In this case it turns out that the physical inclusion of the flux in the manifold matters 
beyond simply shifting the azimuthal angle, and a fractional localized charge might be 
observed. 

This potential cross-over is supported by the results of reference [31 J which 
numerically solves for the case of a constant magnetic field through the two disjoint 
surfaces: a cube with x and y faces identified and open z faces. They find the expected 
fractional charge smoothed out over each face, is preserved. 
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Figure 10. (Color online) Cross-section of a configuration supporting a qubit. The 
blue-green band surrounding the topological insulator represents a ferromagnetic mass. 
See text for detailed explanation. 



5. Flux-Charge Qubit 

The pattern of the results from the previous section implies that unit fluxes induce 
new states with energy equal to the local mass (as opposed to the Hall current which 
depends on the masses globally). Taking this into account, it is possible to conceive 
of a configuration that will produce close to degenerate E — states separated by 
an arbitrary energy gap-a desirable property for a qubit. Qubits using topological 
insulators have been proposed in reference [32] involving Majorana fermions. The 
current proposal would use (fractional) electric charge to distinguish between the states, 
while requiring magnetic flux for their stabilization. It might therefore be considered a 
'flux-charge' qubit. The basic idea is shown in figure [lOj with details described below. 
One imagines the cylinder coated with mass almost entirely (shown in blue around the 
topological insulator surface) except for a region of radius r 2 . If r 2 ~ 1/m then any state 
localized in r 2 will be gapped by order m as well, thereby making the entire spectrum 
gapped by m. Now if a localized flux can be inserted within r 2 , with width of say ri, it 
will induce a low energy states which in the limit of an ultra-localized flux will approach 



E = (two such states localized on the top and bottom windows as shown in figure 12). 

While in previous sections I have been considering a single unit of flux, when more 
than one unit of flux is inserted, other states which start at lower k flow down into the 
gap towards E=0. The flow is along the same, or similar, trajectory. A generic picture 



is that of figure [TT] which shows the flow of energy states as up to two units of flux are 
inserted, after which there will be a new pair of states with E ~ ±0 at k = 3/2, and 
k = 1/2, while the rest of the spectrum remains at \E\ > ~ m. Each of these states 
will be localized in different regions: E ~ +0 states will be localized in the bottom 
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Figure 11. (Color Online) The energies of the induced state as a localized flux 
is inserted for the configuration of figure 10 The energy flow for up to two units 



of flux are shown. The four curves from top to bottom are for different flux radii: 
7*1 = 0.5r2, 0.1r2, 0.05r2, 0.001r2 respectively. All cases have T2 = 0.1R and m = 10 
R _1 . States starting at k = —1/2 flow to k = 3/2, while a second state starting from 
k = —3/2 flows to k = 1/2 along similar curves. Note that symmetric E < states 
and states that have energies \E\ ~ 10 (as in figure 12) are not shown. See section [6] 
for discussion on magnitudes. 
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Figure 12. (Color Online) Electron occupation (filled points) and insets showing 
charge density (plotted on the surfaces) for qubit states |1) (left) and |2) (right). 



Solutions are for the same parameters as those in figure 11 



flux-piercing while E ~ — states will be on top. Therefore E > and E < states 
will not be mixed by local noise (different k states however can be). If, however, the 
flux is fractional there will in general be one state in the mid-gap region spoiling the 
energy separation (between k = —1/2 and 1/2 in figure 11). Therefore a limiting factor 
is the need for integer flux. A second issue, is that for the energies to flow to a value 
sufficiently close to E = 0, ri, the flux width, must be sufficiently small compared to r 2 , 
otherwise the mass exerts a local proximity effect for the k = 1/2 state. From figure 



11 



at least r\ < 0.1r 2 is required for reasonable energy scale separation. 

With interactions turned on, the true ground state is the most neutral occupation, 
which is zero for even flux and ±e/2 for odd (split on caps). A fluctuation through an 
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even flux would destroy the charge and so for all practical purposes precision up to a 
single flux is needed, or specifically exactly one flux quantum. 

These constraints may have a potential solution. To shield flux from regions where 
it is not wanted, a (type-I) superconductor can be put over the top and bottom 
surfaces as shown (figure 10) with holes of radius r\. The induced current of the 
superconductor would favor integer flux through ri, which because it is small could 
easily be accommodated by a reasonable macroscopic magnetic field over the whole cap. 
To further stabilize the single flux constraint, a Type-II superconductor could be used 
with the temperature of the system tuned so that the vortex coherence length ~ r\ 
favoring a single vortex. Other possibilities perhaps using SQUIDs might also be suited. 
If a single flux quantum can be stabilized then the basis of the qubit are simply 



the occupation of E ~ 0^ + ^ state, as shown in figure 



|2> 



l 0+ 



Y\.E<m 



12 



1 



n 



E<m 



|o>, 



Appendix A). |1) has a — e/2 charge at 



|0) (in the convention of 
the top window (+e/2 at bottom) and vice versa for |2). Then in this basis, applying 
an electric field, or a potential V with V\ top = —Vbottom = v, generates a term in the 
Hamiltonian ~ va 3 in the qubit subspace. A a 1 or a 2 matrix element necessary for 
full unitary evolution, is less trivial. Effectively switching the electron occupation is 
required. One potential route is simply inserting metallic electrodes connecting the top 
and bottom windows, allowing the electron to tunnel and with the aid of a bias voltage. 
Exploring these possibilities are left for future work. 



6. Discussion 



In summary the details of how anomalous currents manifest themselves on the surface of 
topological insulators has been presented as well as a potential application in the form of 
a flux-charge qubit. Most striking is the extended proximity effect, which is nonetheless 
rendered reasonable by a lack of quantization. This effect is manifested through the 
non-local nature of the wavef unctions, despite the spectrum responding to mass terms 
locally. The lack of quantization is more physical than the sharp Sign(m) dependence 
found for the infinite plane. The existence of a chiral band which was previously known, 
is consistent with the result that oppositely orientated Hall phases are fully quantized, 
so that Chiral band receives or releases a whole charge. An effective unquantized current 
is restored in the sense that the chiral band can be diffuse. 

It was noted for the torus that mL/2ii ~ 20 was necessary for the onset of full 
quantization while for the closed cap geometry, from figure [8j onset occurs at mR ~ 1/8 
for massive caps and md ~ 8 for mass on the sides of the cylinder (case c). For a typical 
Fermi velocity fT4l [T5l [T6] Vf = 10 6 m/s, setting L, rf, or R on order of 1 mm would 
correspond to m = 8 x 10 -2 , 5 x 10 -3 , 8 x 10 -5 meV respectively, or for L, d ) R of order 
1 /im, full quantization begins at m — 80, 5 and 8 x 10 -2 meV. These are well within 
experimentally observed time-reversal breaking surface mass terms of order meV [19J. 
The wide range of values suggests that the effect is strongly geometry dependent, and 
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as such these values should only be taken as rough guides. 

The required magnitudes for the qubit proposal are more constraining. Using the 
value for Vf above, for a gap of 3° k ^ 0.3 meV, a value of r 2 = 2.5 fim is required. 
Therefore r\ < 0.1 /im, which is small. The value of R and the aspect ratio is essentially 
unconstrained. 

Finally, whether fractional charges remains intact for a diffuse enough magnetic 
field remains an open question which can be investigated both numerically and 
experimentally. If this is the case, then the flux-charge qubit could be supported and 
offer an alternative approach to quantum computation. 
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Appendix A. Spectral asymmetry and anomaly 

The quantum anomaly can be understood by two main routes. For a (2+1) dimensional 
Dirac theory on the infinite plane coupled to a background gauge field A^{x) ) one can 
simply compute the effective action and extract the current term j2j [I] 



S 

(jn) = i—— In det(i0 + eA - m) (A.l) 

= ^rV^ + --- ( A - 2 ) 

8tt \m\ 

••• (Q) = J <J°> = \* (A.3) 

where the field strength F^ y = d^A y — d u A jjL and $ = — e J d 2 rB / {2ti). While this gives 
an exact analytic result, explicit computation of the log of the determinant is difficult 
in all but the simplest cases such as the infinite plane. 

Alternatively one can diagonalize the single particle Dirac Hamiltonian H ) 
Hipx(x) = E\ipx(x). are n-component classical spinors, H contains 7 M matrices 
satisfying {7 M (x), 7^(x)} = g^(x) and g is the coordinate metric (see references [331134] ). 
Then the fermion annihilation operator at x can be expanded in this basis [35j : 

) = a °^° + Y a +A^+A(x) + b -\*p-\(x) ( A - 4 ) 

o +A —A 

where ip±\ are the positive (negative) energy eigenstates, ip Q are possible zero energy 
states, and a, b are fermionic annihilation operators satisfying the usual commutation 
relations (here and below, A with no ± qualifier is taken to mean any eigenstate). 

Then the vacuum expectation of the current, J^(x) = ^\^\x)rf^^(x) — 
(^(x)) Ts (^' f (x)7 7 M ) T5 ], (with 7 7(2?) and i/j if)(x)) can be written [[35]: 

±A 
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+ 1 E V>hVV> -\Y1 V>hVifc. (A.5) 

O It 

where ±A refers to non-zero states and o, u sum over occupied(unoccupied) zero modes. 
To extract the charge in one region versus another (J°(x)) may be integrated over 
regions of interest. 

Appendix B. Wavefunctions and quantization condition on the torus 

Although it is simple enough to write the Dirac equation and the form of solutions for 
both of the cases in the text, for clarity in notation, I will split them. For a review on 
vielbeins and spinors in curved space see references (36J [371 [34] . 

Appendix B.l. Case 1: Fully massive torus; step potential in (ft direction 

As discussed in the text (with unit definitions), for a fully massive torus with a potential 
V a = (— v, +v, —v) for regions I, II, III as in figure [l] the Dirac equation is: 

{-io x d z - io 1 d 4) + a 3 m)ip a = (E- V a )^ a (B.l) 

To simplify the formulae I have put v instead of — ev/2. In this case the z direction is 
trivial and solutions are of the form ijj = exp(ikz)f a ((j)) with k e Z (the same k in all 
regions). If E ^ ±\/k 2 + m 2 ± v and ^ ±|fc|, / a (0) is given by: 

f(fh\ - A ( , 1 I P W(E-V a y-k^-m^ 

J VP) - A* I iy/(E-V a ) 2 -k 2 -m 2 +k J e V 

V E-V a +m / 

-+- b I , — o i e -W( E - v ^ 2 - k2 - m ^ (R 2^ 

+ n ot I -i^/(E-V a ) 2 -k 2 -m 2 +k ] e V \ n - Z ) 
\ E-V a +m / 

The coefficients in each region and energy quantization are determined by normalization 
and matching conditions fi(L 2 /2) = V f m (-L 2 /2) : h{l/2) = f u (l/2) and /m(-//2) = 
fn(—l/2) with r] = — 1 for the chosen spin-structure in the cj) direction. The quantization 
condition is: 

|a| 2 (cos(p/ + q(L 2 - I)) - 1) = |6| 2 (cos(p/ - q(L 2 - I)) - 1) (B.3) 
where I have defined q + v) 2 — k 2 — m 2 , p — v) 2 — k 2 — m 2 , a 

E-vZrn - and b = E-vlm ~ This is solved numerically and shown (for a 



representative case) in figure [2J The special cases E ^ ±^k 2 + m 2 ± v and E ^ ±|fc| 
(and special sub-cases of these such as E = 0) must be treated separately although in a 
similar manner as above. The solutions in those cases are linear in one of the regions, 
and in general do not yield new solutions. 

Appendix B.2. Case 2: Partially massive strip 
In this case the Dirac equation is: 

(-ia 1 d z - ia 2 d^ + a 3 m a )ijj a = E^ a (B.4) 
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where m a = (0, m, 0) for regions I, II, III as in figure g Different spin-structures are 
encoded in whether the boundary condition on ip is anti-periodic or periodic in (j). The 
In each region I, II, III the solutions are of the form ip = exp(ik(j)) f \z) with k G Z + 1/2. 
If E 7^ ±\/k 2 + m 2 and E ^ ±|fc|, /(^) is given by: 

/(*) = 




(B.5) 

E+m / 

where the region label has been omitted and the mass m is understood to be zero in 
region I and III. The coefficients in each region and energy quantization are determined 
by normalization and matching conditions /i(L 2 /2) = r]fiu(— L 2 /2) : f\{l/2) = fu(l/2) 
and /m(— 1/2) = fu{—l/2) with rj = 1 for the chosen spin-structure along z loops. The 
quantization condition is as before: 

\a\ 2 cos(pl + q(L 2 - I)) = \b\ 2 cos(pl - q(L 2 - I)) (B.6) 

but in this case the definitions are q = \/ E 2 — k 2 ) p = V ' E 2 — k 2 — m 2 , a = — , 

and b = — This is solved numerically and shown (for a representative case) 

in figure 5, The special cases E = ±^k 2 + m 2 and E — ±|fe| (and special sub-cases of 
these such as £" = 0) must be treated separately although in a similar manner as above. 
The solutions in those cases in general do not yield new solutions. 

Appendix C. Derivation of the Dirac equation on closed cylinder surface 

In the notation of figure [6| the closed cylinder is divided into three regions. The local 
form on the caps is as the Dirac equation for the disk described in equation [6] which 
can be gotten by transforming the Dirac equation in Cartesian coordinates into polar 
coordinates. To guess the form on side, region II, one cannot simply replace r — >> R as the 
metric suggests. In particular the spin connection is wrong, and stems from the fact that 
all the curvature is at the sharp corner where the coordinate transformation is required. 
For a review on vielbeins and spinors in curved space see references [36, EH E3]. A real 
manifold must be smooth and the transition between charts should be a map between 
IR 2 IR 2 , although in practice one dimension can be made of infinitesimal width if the 
manifold is smooth. This suggest a more careful procedure is to smooth the corners into 



a semi-circle of radius 5 as in figure CI. Replacing the relation 5 2 (r — R) 2 + (z + dy2 



2; 



into the 3-dimensional Cartesian metric, the induced metric is, 

S 2 

9 — } ^— rrfr ®dr- r 2 d(f) ® d(f) (C.l) 

o z — (r — R) z 

and leads to the Dirac equation (using (j) dependent vielbeins) in (r, (p) : 

-^1 . k ,VQ3 4 _ k a (C . 2) 
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z = -d/2 



r=R 

Figure CI. (Color online) Side view of smoothed out cylinder edge. 



2or 



(C.3) 



As r i?, the equation becomes 

just what we expect on the cap evaluated at r = R. However, replacing r(z) and taking 
z — >■ =n one gets: 



(C.4) 



e 1 2 



ia x d 7 



dch 



(C.5) 



+ 

just what one expects from the side using an embedding, without any further 
transformation of ip i. e. i/Jm(R,(f)) = < 0n(— d/2, 0). Less trivially is the boundary 
condition between the top and the side. The derivation proceeds in the same way 
if the local z were facing down. So if one defines a new z n = —z, and calls 

■a 3 

the tentative wavefunction ipu n with the notation ipu n = e^^nn one would have 
found (— io x d Zn — i(J 2 ^)^ii n {z n) (j)) = E^u n (z ni (/)) (with boundary condition ^i\ r =R = 
^iin\z=d/2) which is (+ia 1 d z — icr 2 jj-)ijjii n (z,(f)) = E^u n (z^). Evidently, then, the new 



ihin — v ipih since a a a 



In summary equations M and 11 are obtained. 



Appendix D. Wavefunctions and quantization condition on cylinder 

The local solutions to equation 9^ (if there is no flux) are for the caps ^nii 
exp^WWr). If E , ± J^T± and E + ±m , m , /M ha, the form: 



J 4 As/W 



m*r) 



m-c 1 \ l ^ Lj vw — ?r) 



G 



Y~ k _i_(y/E 2 -m 2 r) 
2 1 *V^V (V £23^2 r) 



(D.l) 
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where a region label =1, III should be understood for the mass m, f(r) and the 
coefficients C\ and C2. /cEZ + l/2is globally the same value for all regions, and 
different k are linearly independent as required from the boundary condition along the 
(j) direction (they are good quantum numbers). J n (r), Y n (r) are Bessel functions, the 
two independent solutions satisfying r 2 ^ + (r 2 — n 2 )h — 0, J, Y — h(r) ( J n (r) and 
J- n (r) also work if n is not an integer). The solution on the side is ipu = exp(ifc0)/n(r), 



and if £" 7^ ±y k 2 + and E 7^ ±mn, /(r) is (again omitting region II label): 

/(*) = 




_^ ; (D.2) 

Square-integr ability at the origin sets one of the coefficients of the cap to zero, 
when |fc| > 1/2. Using the boundary conditions discussed in the previous section, the 
remaining coefficients are determined along with the quantization condition: 

-2^ 2 -^-m? ; ( ihiig* - gm } = ign* ~ hi ^ ^ 



gin + ih n n J hi + ign 



with h a = ^^i(^ 2 -m^), g a = J~ k _ h (^E 2 -miR) and 7 

^E 2 -k 2 -m 2 1 +ik 



E+mu 



2 v - ,v 2 

When a localized flux is inserted through the origin, then away 



from origin, the solutions are still of the form given by (D.l) with the replacement 



k k + $ for $ flux quanta. In this case, square integrability at the origin 



no longer constrains the coefficients in (D.l), instead a matching to the solutions 



in the B-field region must be done. As discussed in the text, the simplest flux 
profile allowing for analytical solution is to take a delta-function ring [29]: B{r) = 
<M(r — e)/(27re). Then inside the ring, normalizability again constraints one coefficient 
to zero. Matching this at r = e gives a required relation between the coefficients 
in the interior of the flux and the rest of the cap. Then the quantization condition 



is approximately the same as (D.3) with k —> fc + $for \k + $| > 1/2. For 



\k + $| < 1/2 the explicit form of the B-field must be used to match wavefunctions 



at e. The result can be summarized by making redefinitions in equation D.3| (assuming 



* > °) : K = ^^(d 1 J kH+ ^E 2 -m 2 a R) - d 2 J_- k _ h ^E 2 - g a = 

diJ~ k _i +q> (^E 2 - m 2 a R) + d 2 J_ k+ i_< s> (^/E 2 - m 2 a R) where 

d 1 = J i-w- d2j --w-* d 2 = J £-i/2 + «Jfc + i/2-Ji + i/2 + *J£-i/2 with all unspecified 

J fc-l/2+<£ J -fc-l/2-3> J fc-l/2+<£> J -fc + l/2-3> J fc + l/2+<£ 

Js evaluated at \/E 2 — M 2 e so long as $ is not an integer. For $ integer, one further 
must replace J _~ k to Y (the second Bessel function) and — d 2 to +rf 2 in h a . 
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